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We study the basin of attraction of static extremal black holes, in the concrete setting of 
the STU model. By finding a connection to a decoupled Toda-like system and solving it 
exactly, we find a simple way to characterize the attraction basin via competing behaviors 
of certain parameters. The boundaries of attraction arise in the various limits where these 
parameters degenerate to zero. We find that these boundaries are generalizations of the 

(-H 

recently introduced (extremal) subtracted geometry: the warp factors still exhibit asymptotic 
integer power law behaviors, but the powers can be different from one. As we cross over 
one of these boundaries ("generalized subttractors"), the solutions turn unstable and start 
blowing up at finite radius and lose their asymptotic region. Our results are fully analytic, 
but we also solve a simpler theory where the attraction basin is lower dimensional and easy 
CN to visualize, and present a simple picture that illustrates many of the basic ideas. 
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1 Introduction 



In this paper, we will consider attractor black holes [H EJ [31 H] : specifically, we will study 
the basins of attraction of static non-supersymmetric attractors mostly in STU supergravity. 
Attractor black holes are extremal black holes that can exhibit scalar hair for a fixed horizon 
value of the charges. Since the thermodynamic properties of the black hole are determined 
by its charges, the attraction phenomeon is a hint that only the near-horizon geometry of 
the black hole might matter for understanding its thermodynamics. 

In a recent paper [5] , we showed that the recently introduced subtracted geometry [6j [7J [8] 
(in its extremal limit) can be thought of as a boundary of an attraction basin. Subtracted 
geometries are black holes in Einstein-Maxwell-dilaton theories that arise upon replacing a 
certain warp factor in the metrics of flat space black hole solutions by another (specific) warp 
factor. The precise algorithm for the replacement depends on the original black hole metric, 
and the replacement changes the asymptotics. But it does not change the black hole's 
thermodynamics and it manifests a certain "hidden" conformal symmetry of the original 
black hole. The original motivations for introducing the subtracted geometry can be found 
eg., in [9} [Hll [HI [El [131 El [151 HE! [T7j, but we will not concern ourselves with that here. 
Our goal will be to study the attraction basins of some exactly solvable attractor systems 
analytically, and to identify the various boundary behaviors that emerge. We find that the 
extremal subtracted ( "subttractor" ) geometry that was found in [5] and some generalizations 
of it arise naturally as boundaries of the attraction basin: the original subtracted geometry 
of j6j [7J [8] had a warp factor that went linearly in r for large r, but we find that in the 
STU model, attractor boundaries can have behaviors that go as ~ r, r 2 or r 3 . We call them 
generalized subttractors. They can all be understood in terms of replacements of the warp 
factor in the original geometry. 

The systems we consider are fully integrable in terms of a Toda system, so we can 
find exact solutions and get a complete understanding of the attraction basin. We can 
characterize the general solutions in terms of certain integration constants d a . The solutions 
are regular everywhere at and outside the horizon only when d a > 0. Therefore, the entire 
attraction basin is captured in <i a -space by the first quadrant for a simple one-dilaton theory 
where d a can be d\ or c? 2 and by the first orthant in the STU model where a can take four 
values. It will be interesting to see what generalizations of these statements can be made for 
other Einstein-Maxwell-dilaton theories, but the question is not necessarily straightforward: 
we got lucky here because the system is Toda integrable. 

We find that the boundaries of the attraction basin occur when one (or more) of the 
daS degenerate to zero. We provide a full characterization of such boundaries. The unstable 
solutions corresponding to negative d a do not have an asymptotic region and diverge at finite 
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radius. 

The plan of this paper is as follows. In the next section, we describe our system: the black 
holes we consider all fall into a certain attractor ansatz (even though all the black holes al- 
lowed by the ansatz are not extremal/attractive). We show that under certain circumstances, 
the coupled system of Einstein-Maxwell-dilaton equations of motion can be brought to the 
form of a generalized Toda system and can be integrated exactly. The systems we consider 
fall into this category, and can be diagonalized and integrated. In section [3j we introduce 
our main object of interest, namely static black holes in a consistent truncation of STU 
supergravity. We solve the system completely using our Toda approach and describe the 
attraction basin and the boundaries in detail. Since the STU system has three scalars and 
the attraction basin in terms of the d a parameters is four dimensional, we then turn to a 
simpler system. This system has the advantage that some basic ideas remain the same, 
yet we can get a more intuitive feel for the attraction basin. The theory we consider can 
be obtained from the STU system [5] and is also integrable and we find a fully consistent 
picture. 

The solution space is quite rich and we find various auxiliary results on the way but we 
have relegated them to appendices. Among these are the observation that the "vertex of 
attraction" is an AdS2 x S 2 space, and the identification of the slice along which a given 
black hole flows in the attraction basin when perturbed. We conclude with some comments 
about various connections and open questions. 



2 A General Toda-like System 

A general Lagrangian arising in (ungauged) supergravity is 

S = j d A xV=~g(R - 2m) 2 - f ab (4>i)F^ u F b - \u(<M)F^F%^) (2.1) 

where % = l,...,n is the number of scalars and a,b = 1,...,N counts the number of U(l) 
gauge fields. The /„& and f ab are the gauge couplings including the axionic piece. This is the 
form considered in PQ. We will look for solutions of this system in the "attractor ansatz", 
which means we will consider solutions of the form 

dv 2 

ds 2 = -airfdt 2 + — -r + b(r)W, (2.2) 
a(r) z 

F a = Q a m sm6d6 A d<t> + L ^ l (Q eb - f hc Q c m )dt A dr, & = &(r). (2.3) 

x We emphasize that the solution has previously been found in a different form in |18j . but we will describe 
some advantages of our approach when the context arises. 
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The equations of motion take the form 



{a 2 b 2 )" -2 = (2.4) 
b" 9 

j + = (2-5) 
and the first order constraint from Einstein equations is 

2i/2 , 2'j,2' , ^eff(^t) 2i,2 jl/ 2 in /<, 7 \ 

a o + a o H — a o q>i —1 = 0. (2.7) 

cr 

The effective potential that shows up in the equations of motion is 

VMi) = f ab (<f>i)(Q ea - hMi)Q c m )(Qe b - hMi)Q d m ) + f a MQ a m Q b m - (2.8) 

At this stage, generalizing pQ we introduce 

Ui = (pi for 2 = 1, ..,n, u n+ i = logo, z = logab. (2.9) 

With these definitions, via simple linear combinations, one can bring the equations to the 
form 

u n+ i = e 2 ^V eS { Ui ), (2.10) 

Ui = \e 2u ^d u y cS ( Ui ), (2.11) 

z = e 2z . (2.12) 
The energy constraint is redundant and takes the form 

u 2 + u 2 n+1 + e 2z -z 2 = e 2u ^V eS (u t ). (2.13) 

The dots stand for derivatives with respect to the tortoise coordinate 

/dr 1 / 2m\ 

^=2^ l0g ( 1 -v) < 2 ' 14 > 



T 



We choose boundary conditions such that the solutions are regular at r = 2m. This is what 
we call the outer horizon and we choose our notations so that they match with the black holes 
in [SI E3] in the appropriate limits. The inner horizon has been fixed by a coordinate choice 
to be at r = 0. The energy constraint has been used in bringing the first three equations to 
the simple form presented here. The equation for z is really the statement that (a 2 b 2 )" = 2 
and is solved by 

b 2 = r(r - 2m)/a 2 , (2.15) 
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once we make our coordinate choices. So (2.10 2.11) are what we really need to solve 



So far our set up is completely general within the attractor ansatz. To proceed further, 
we will assume a less general form for V e s- The form we will consider will still be general 
enough to contain the previous exact solutions [TJ |5] and as we will show in the next section, 
the STU attractor also falls in this class: 



N 

,i\MI,) '- } J 

a=l 



V P fe(Ui) 



(2.16) 



As explicitly indicated, there is one summation over a and another one over i. Now, we will 
assume that 



N — n + 1. 



(2.17) 



The notation in the following will be that summation over repeated indices (irrespective of 
how many times they are repeated or the precise index placement) will run from 1, N for 
a, 6, c, d and from 1, n for i,j, k. Using this, and defining^] 



a a b = (a a )b = (Olai «a7v) = (otai 2^) 



(2.18) 



brings (2.10 2.11) to the form 



U; 



2 pOiabUb 



6«6 



Ql e 



(2.19) 
(2.20) 



With our choice of the effective potential and because of the relation (2.17), these two 
equations can be combined to the form 



u 



- A D 2 p a cb u b 



(2.21) 



where the square matrix A is defined via 



(2.22) 



and a is the number of the column. Defining the matrix m a b 
change of variables 



X a = (A l )abU b +{m l )ab^Q 



a ac A c i, and making the 



(2.23) 



2 We can (and will) think of a ab as a square matrix a or as a set of row vectors a a depending on the 



context. In (2.18), 2 a stands for a row vector indexed by a in which each slot contains the number 2. 
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brings us to the final form 

X a = e mabXb . (2.24) 

This is a set of N differential equations which take a simple structure: they are coupled to 
each other linearly in the exponent. The matrix is a constant matrix, and therefore our 
equations can typically be related to Toda systems via simple tricks (see eg. [U El [191 120])- 
We will do this for static solutions in the STU model in the next section. 



3 Toda Solutions in the STU Model 



We will mostly be working with a theory that has three scalars and four vectors in most 
of this paper, that arises as a consistent truncation of TV = 2 supergravity with four vector 
multiplets El 1230 



J d'x^-g[R- l -{d Vl f-\{t 



(3.1) 



We look for an attractor ansatz for this action of the form presented in the previous section, 
with the specific choice for the gauge fields: 

Q 2 



F 1 = 2Qx sin 9d9 A 



F- i = 7 1 = 2Q 3 sin 9d9 A 



4g2(01-</ , 2-03) 



7 l 



2b(r) 2 

4e 2(^. 1 +^ 2 +0 3 ) 



dt A dr, 

dt A dr, 



Qa 



2b(r) 



(3.2) 



(3.3) 



This is of the ansatz presented in (|2.3|). With these, the equations of motion take the form 
the ] 

ell = Qi< 



presented in the previous section, with ipi = 20 i; and with the effective potential 

y pff = Q2 e -20 1 +2</, 2 -2</. 3 _|_ g2 e 2</. 1 -20 2 -2</» 3 + g2 e -20 1 -2 < /. 2 +2^ 3 + g2 e 2</. 1 +20 2 +2</. 3 ^A) 



This structure falls in the general form for the effective potential that we presented in the 
previous section. In fact with the choice 



«i = 01, «2 



«3 = 03, u A = \oga 



(3.5) 



the matrix 



A 



( 


-1 


+1 


-1 +1\ 




+1 


-1 


-1 +1 




-1 


-1 


+1 +1 


\ 


+1 


+1 


+1 +1 / 



(3.6) 



3 We can also think of the system as arising in the STU model, up to an electric-magnetic duality. At the 
level of the effective attractor ansatz theory that we write down, this makes no difference. 
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and the a a b matrix can be read off from it. This results in the m a b matrices becoming 
diagonal, and therefore the system decoupling into four separate ODEs: 



X a = e 2X ». (3.7) 



This is a trivial version of the Toda lattice equations where there is no coupling between 

dv 

dx- 



neighbors and can be immediately solved by noting that X = v 4^, where v = X. The end 



result is 

^(ssr^)^(l)- ^ 

where c a and d a are integration constants (no summation over repeated indices here) and we 
have defined the F a in terms of the sinh functions for future convenience. This is the most 
general solution of the attractor ansatz in the STU model with effective potential given by 



The above decoupling was observed in a different, but equivalent set up in the STU model 
in [18]. It seems they came to this result by inspection and cleverness, we see here that it 
can be accomplished by cranking our Toda-ization machinery. The reason why we prefer 
our way here is that the integration constants that naturally arise in our set up have an 
immediate and clean interpretation in terms of the attraction basin. Of course, even though 
the forms of our solutions look superficially different, we show that the [18J solution can be 
brought to our form by connecting our parameters to theirs carefully ( ]3.20||3.22 ). 



Translating back from the Toda-style variable X a to the original variables, we find that 

Q2Q4 c lC3 F 2 F 4 ' 1 ' ] 

A<j> 2 _ ^2 Q3 c\ c 4 F 2 F 3 

' ' Qi Qi c 2 c 3 F 1 F 4 ' [6 - W) 

e 4<fe = Q± C 3 C 4 ^1 F 2 ^ ^ 3 U ^ 



■2 I ci c 2 c 3 c 4 1 2 r(r - 2m) 



a V2 4 QiQ 2 Q3Q4F 1 F 2 F 3 F 4 ' b a 2 • (3 ' 12) 

These solutions remain invariant under the c a -H- —c a , so we will work with c a > 0. We 
want regularity as r — > 2m (or, r — > —00). The finiteness of the scalars yield 

ci + c 3 - c 2 - c 4 = c 2 + c 3 - c\ - c 4 = c\ + c 2 - c 3 - c 4 = (3.13) 

which forces 

ci = c 2 = c 3 = c 4 = c. (3.14) 
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The finiteness of b at the horizon (which is necessary because otherwise the horizon will be 
singular) then gives us 



m. 



So the solutions which are regular at the horizon are 



,403 



Qi 


Q 3 


sinhm 


V 


-di) 


sinhm 


V 


-da) 


Q2 


Qa 


sinhm 


V 


-d 2 ) 


sinhm 


V 


- d A ) 


Q2 


Q 3 


sinhm 


V 


-d 2 ) 


sinhm 


V 


-d 3 ) 


Qi 


Qa 


sinhm 


V 


-di) 


sinhm 


V 


- d 4 ) 


Qi 


Q2 


sinhm 


V 


-di) 


sinhm 


V 


-d 2 ) 


Q 3 


Qa 


sinhm 


V 


-d 3 ) 


sinhm 


V 


- d A ) 



m 



(Q1Q2Q3Q4 



-1/2 



(3.15) 

(3.16) 
(3.17) 
(3.18) 

.., . . . _ : . , (3.19) 

2 w sinh m(r — di) sinh m(r — d 2 ) sinh m(r — d 3 ) sinh m(r — rf 4 ) 

These solutions are general. In particular, we haven't imposed asymptotic flatness or ex- 
tremality at the horizon. We can show that when all the <i a 's are positive, the solutions are 
regular everywhere outside the event horizon (which also is regular), and asymptotically flat. 
But first we compare this to the solutions of [18] . Remembering the relation (2.14) between 
r and r, it can be shown that our solutions turn into theirs (see equations (2.32-2.35) in 
[IS] ) when we set our scalars (pi = T]i/2 (where r]i are their scalars) and choose 

Qx = £3/2, Q2 = £i/2, Q 3 = B 2 /2,Q 4 = q /2. 
2 sinh md\ 2 2 sinh md 2 2 



cosh mdi — sinh md\ 
2 sinh mdz 



a 



"' coshmc/2 — sinh m<i 2 
2 sinh mdi 



(3.20) 
(3.21) 



a 



cosh mdz — sinh md 3 ^ cosh mdi — sinh mdi ° 
where the right hand sides are in the notations of [IS] and the left hand sides are our 
notations. The superficially strange assignment in the index matching arises because the 
theory has a cyclic symmetry in the a = 1, 2, 3 (but not 4) index and their choice is related 
to ours up to a cyclic shift. Another comment is that their gauge fields are different from 
ours in the choice of the electric-magnetic duality frame, but the effective attractor theory 
is the same either way. 

For our purposes, our form of the integration constants is much more convenient. This is 
because as we will show presently, they are directly related to the data at asymptotic infinity. 
The [IS] integration constants involve m which is an IR (near-horizon) piece of data. This 
mixed UV-IR notation for the data makes things less convenient to understand and work 
withn so we will stick to our Toda-notation. 



4 in particular, note that the extremal limit involves m — > and the coefficients do, ■■■0,4 in |18j collapse 
to zero in that limit. This indicates that these are not a good notation for working with attractors. 
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Now we will demand that the geometry is asymptotically flat and that the scalars take 
finite values at infinity, ie., we demand 



a — ► a , <pi — > <p ioo (3.23) 

We can use these to fix d a easily enough in the non-extremal case, but since we are interested 
in attractors, we will do this in the limit m — > 0. We get three equations from the scalars 
and one from a(r), so we can fix the four integration constants via 

di = e + ^ °-^ ° +</,3o °, (3.24) 



2 a \Q 
1 

2 a |Q 2 | 
1 



2 a \Q 



3 



g-</>loo+02oo+^3oo 25) 
g+01oo+^2oc— </>3oo 26) 



d A = - — l —— e-^- (3.27) 
2 ao Iv4| 

The absolute value signs are there to ensure that if any of the Q a 's are negative^] the corre- 
sponding d a will have an overall negative sign in its expression. This is necessary to make 
sure that all the c? a 's are positive, which is a necessary condition we will further elaborate 
on in the next section. 

We have left the boundary value of ao to be an arbitrary (positive) constant at this 
stage, but in the attractor system every term that contains a also contains a derivative of 
r, so the asymptotic value of the metric can be rescaled to unity by absorbing it in r. So 
these are asymptotically flat solutions. The fact that the asymptotic value of the scalar 
can be arbitrary is a direct consequence of the attraction phenomenon. Now we turn to a 
characterization of the attraction basin. 



4 Basin of Attraction 

Even though the asymptotic values of the scalar can be arbitrary, the scalar perturbations 
at the horizon cannot be arbitrarily large. When they become too large, the solution stops 



flowing to asymptotically flat solutions, see Appendix D Indeed, we will show that such 
solutions do not have an asymptotic region. This is what characterizes the boundaries of 
the attraction basin. We will try to understand them in this section. 



5 Note that the possibilities for the signs of Q a 's are such that one can have all of them positive, all of 
them negative or any two of them negative while the other two positive. Other choices are not allowed. 
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In the extremal limit (m — > 0), our solutions take the form 

(4.1) 
(4.2) 
(4.3) 

V3 V4 1-l T u, 3 i ) yi t 114:1 ) 

a2= i {QiQiQzQi)- 112 r 2 

" 4 ^(1 + d ir ) (1 + d 2 r) (1 + d 3 r) (1 + d 4 r) ' 

A crucial observation is that if we demand that a is regular for all values of r G (0, 00), then 
we are forced to choose all d a > 0. For negative values of any of the d a , there will be a 
divergence atr = — l/d a > 0. 

The sign choices of Q a 's result merely in some trivial changes in the solution, so we 
will assume in what follows that all the Q a 's are positive. In particular, note that the 
non-negativity of the daS is independent of the signs of the QaS. 

We claim that the boundaries of the attraction basin correspond to some of the d a 's 



e Hi = 


Qi 


Q 3 (1 + d ir ) (1 + d 3 r) 




Q2 


g 4 (l + d 2 r)(l + d A r) ' 


e ^2 = 


Q 2 


Q 3 (l + d 2 r) (l + d 3 r) 




Qi 


Q 4 (1 + dir) (l + d 4 r)' 




Qi 


Q 2 (1 + d x r) (1 + d 2 r) 






Q A {l + d 3 r) (l + d 4 r)' 


1 




(QiQ2QM- 1/2 r 2 



degenerating to zero. It is clear from the forms of the scalars in (3.24 3.27), or from taking 



the r — > 00 limit of the (|4.l||4.4|) equations, that these correspond to the boundary values of 
the scalars diverging^ 



We classify the various possibilities below. 
<) c?4 7^ case. 

• di = d 2 = d 3 = corresponds to <pioo, 02oo, 03oo — > —00. 

• d\ = d 2 = 0, d 3 7^ corresponds to 0ioo,02oo fixed (or diverging slower than 
|03oo|), with cj) 3oc —7- —00 faster than the others. The other cases where two of the 
0ioo(^ = 1 ; 2, 3) are zero while the third one is not, can be obtained by symmetry. 

• d\ = 0, d 2 7^ 7^ d 3 corresponds to 0ioo,03oo — > —00, with (fi 2oo — >■ +00. The 
other cases can be obtained by symmetry. 

• None of d±,d 2 , d 3 are zero. It is not a boundary. 
<) c?4 = case. 

• d\ = d 2 = d 3 = 0. This solution corresponds to the "vertex" of the attraction 
basin, and in fact is an AdS 2 x S 2 geometry. We discuss it in an Appendix. 



6 The generalized subttractor boundaries arise as logarithmic envelope curves along which the scalar 
diverges, so there is no contradiction with the fact that the attraction basin has boundaries. 
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• d\ = d 2 = 0, d 3 7^ corresponds to <pioo, 4>2oo — > oo, with 3oo — >• — oo. The other 
cases can be obtained by symmetry. 

• di = 0, d 2 7^ ^ d 3 corresponds to 0ioo,03oo fixed (or diverging slower than 
|02oo|), with 02oo — > +oo faster than the others. The other cases are obtained by 
symmetry. 

• None of di,d 2 , d 3 are zero, corresponds to 0ioo, 02oo, 03oo —> +oo. 

The restriction of positivity on <i a 's implies that asymptotically flat (attractor) solutions 
fall in the first orthant of the (di, d 2 ,d 3 , di) space. The boundaries of the orthant correspond 
to the attraction boundaries. When one of the d 's becomes less than zero, we end up with 
solutions that diverge at finite radius. We will illustrate these matters using a simpler system 
in the next section. 

We present in an appendix the explicit forms of these various attraction boundaries. We 
find that they can all be understood as new solutions with new warp factors (a la "subtracted 
geometry" (6j [7J E]) on the original Cvetic-Youm black hole, but with more general replace- 
ments of the warp factor than the original subtracted geometry. The asymptotic behavior 
of the warp factor can take different powers of r : it can go as r, r 2 or r 3 . The case when it 
goes as r corresponds to the original subtracted geometry of [7]. Note that the structure of 



our extremal solution (4.1 4.4) allows also the asymptotic behavior r° and r 4 . The former 
corresponds to a unique AdS 2 x S 2 solution (which we discuss in Appendix [Bj . The latter 
represents the usual flat space hairy (and therefore highly degenerate) solutions. The in- 
between powers r, r 2 or r 3 correspond to various attraction boundaries with varying degrees 
of degeneracy corresponding to edges (of various dimensions) of the basin. We discuss these 
general subttractors in Appendix [A| 

We note that the standard Cvetic-Youm black hole's j2Tj [221 E3] extremal limit (see 
Appendix where we quote this solution in its extremal limit) corresponds to choosing d a = 
Another interesting case is when d\ = d 2 = d 3 = d^. This is the hairless extremal black 



2|Qa|' 

hole (which is different from the extremal Cvetic-Youm solution). These are black holes in 
the central positively directed ray in the attraction basin spanned by di, d±. 

The structure of the attraction basin in the STU model is pretty intricate because of the 
number of scalars 0j (= 3) and the number of parameters d a (= number of charges = 4). In 
the next section we will investigate a simpler system with only one scalar and two charges 
where the attraction basin is much easier to visualize. We will see that the general picture 
that emerges is very similar. 
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5 Toy Model 



We will work with the model considered in [5j. The action is of the form (2.1), with 



and (ai,a 2 

) = (2^/3,-2/^3). The axionic coupling matrix f ao is zero. The attractor 
ansatz and the equations of motion goes through as before with the understanding that 
the charges are two in number and both magnetic, and the index % corresponds to a single 
scalar. The effective potential in this case also has a structure that enables Toda-ization as 
well as diagonalization [lj^J This means that we can find exact solutions here as well, and 
the general solutions which are regular at an extremal horizon can be obtained by taking 
the m — > limits of equations (4.1-4.3) in [5J. The result is 

r ^/V3 _ 1 Q2 (1 + d 2 r) 

r 2 1 

a 2 = ; . =, b 2 = r 2 /a 2 , (5.3) 

^^QM^il + d^il + d^f 

where 

V OL\) V a, 2 ) 

is a purely numerical factor. The structure here is obviously very similar to the one we found 
in the last section in the STU model, but here there are only two parameters d a = di,d 2 
and only one scalar. As before, the regular asymptotically flat solutions exist only in the 
d a > regime. Unlike the four dimensional space of STU model however, here the d a y s span 
a two-dimensional space and the attraction basin is the 1st quadrant. So it is easy to plot 
useful pictures of the basin here. 

The boundaries are at d\ = and d 2 = 0. By setting a — > ao, 4>i — > 4>ioc, we find 

9-3/8 ql/8 

di = — > e~^ 3 , d 2 = ^= e^^, (5.5) 

a V^\Qi\ a V^|Q2| 

which implies that d\ — corresponds to the upper boundary where </>oo — > +oo and d 2 = 

corresponds to the lower boundary — > — oo. It is easy to check from (5.2) that this latter 

solution corresponds to conventional subtracted geometry in the extremal limit ("subttrac- 

tor"), and it was arrived at in this way in [3]. The warp factor there goes as ~ r. The d\ = 

boundary however corresponds to 0oo — > +oo and it is a new kind of subtracted geometry 

and the warp factor goes as ~ r 3 in the asymptotic region. 



7 This can be understood because this system has a connection with the STU model [5] which we showed 
to be diagonalizable after Toda-ization in this paper. 
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Figure 1: The attraction basin for our toy model when Qi = 7, Q 2 = 4. The attractor value 
of the scalar is ^ log ( t^qi ) — —0.480177. The upper region where the curves diverge 
at finite radius correspond to d\ < 0, d 2 > 0, the lower region where they diverge is where 
d 2 < 0, di > 0. Inside the attraction basin, it is di,d 2 > and the upper (lower) boundary 
is di = (d 2 = 0). The last non-divergent curves are logarithmic in r as r — > 00, and 
correspond to generalized subttractors. The "hairless" black hole corresponds to d\ = d 2 - 
We have chosen to make this plot with d = ^Jd\ + d\ = 1, and then tuning 6 defined by 
di = dcosO, d 2 = dsin8. The different curves in the figure correspond to different values of 
9. The qualitative features of the plot will not change for other numerical values of d and 
more generally, as long as we are using a monotonic curve in the d\-d 2 plane to slice through 
the attraction basin. 
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We would like to scan the solutions for various asymptotic values to get an intuition 
about the attraction basin. The above equations suggest that this is equivalent to scanning 
(non-negative) d\ and &i- However, since the axes of the 1st quadrant in the d\-d<i plane 
corresponds to the boundaries of attraction, there is a whole (non-negative) d\ ray worth of 
solutions corresponding to the g?2 = boundary (and vice- versa). So to slice the attraction 
basin and see both boundaries in an instructive plot, it is crucial that we move in the d\-d<i 
plane along a curve that cuts both (positive) axeg^J In particular, one will not be able to 
see the other boundary by starting with the c?2 = boundary for some value of d\ and then 
cranking up g?2 arbitrarily (while holding d\ fixed). To make the plot we present, we have 
chosen 

di=dcos9, d2 = dsin9, (5.6) 

with d fixed, and then plotted the scalar as a function of r for various values of 9. 

The above choice of the curve was arbitrary. But starting with a given black hole solution 
fixes the choice of d\ and g?2 (one might call these d® and d®), and perturbations around this 
solution are then determined by the leading order scalar perturbation, which we will call (f>\. 
Since a general solution is fixed once we know di and c?2, we can determine the d\ and d<i of 
the perturbed solution in terms of <pi, d\ and d® by matching it order-by-order with the near- 



horizon series expansion of the general solution (5.2 5.3). Since d\ and c?2 are determined in 



terms of one parameter 1; this picks out a unique curve through the attraction basin with 



initial values d\ and d\. We explain this in detail in Appendix D and demonstrate with the 
example where our starting solution is the subttractor geometry of [5], given by a specific 
choice of d\ and d\. 

We discuss some of the features of the plot in the caption of figure 1, and similar state- 
ments hold also for the more general STU system we considered in the previous sections. 
The structure of the attraction basin is such that as one of the d a y s goes through zero and 
becomes negative, the solutions no longer have an asymptotic region where r — > oo. This 
is because (as we mentioned) the (1 + d a r) causes a divergence at r = — l/d a (> 0) when 
d a is negative, and the solutions literally blow up at finite radius. Note however that the 
extremal near horizon geometry and the attractor value of the scalar are still preserved - so 
the near-horizon region is still attractive to the scalars even though they are blowing up and 
losing their asymptotic region^} 



8 Note also that the hairless black hole corresponds to the c?i = d,2 solution, the central positive ray in the 
g?i-g?2 plane. 

9 Since crossing over the d a = boundary brings out this conflicting, self-destructive behavior from scalars, 
one might call them "dramatic solutions" . But perhaps one should not judge these solutions too harshly. 
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6 Comments 



We will conclude with some comments, open questions, future directions, etc. Some of 
these questions are being investigated [24]. 

We have presented generalized subtracted geometries in the extremal limit. These have a 
straightforward generalization to the non-extremal case, even though we have not emphasized 
this in our paper. However, both these results are limited to the case of static black holes. 
What would be more interesting is if the generalized subtracted geometry also exists in 
the case of black holes with rotation. In [7], the warp factor for spinning black holes was 
chosen by demanding that the wave equation is separable. Together with the demand for 
asymptotic linear behavior, this was enough to essentially uniquely fix the warp factor of 
the subtracted geometry. For higher powers in the asymptotic behavior, it seems likely that 
the separability criterion is not enough to uniquely fix the warp factoip^j However, it will 
be interesting to see how constraining separability is. Together with our knowledge of the 
warp factor in the static limit, and by demanding that the equations of motion be satisfied, 
it might be possible to determine a separable, generalized subtracted geometry for rotating 
black holes as well. A related and very likely possibility is that one might be able to generate 
these new solutions via an appropriate Harrison transformation P, 126^] Finally, there is the 
question of constructing attraction boundaries for rotating extremal black holes [28] which 
is a natural extension of these ideas. 

The solutions we have investigated all were integrable via the Toda trick. Furthermore, 
the toy model we considered is secretly connected to the STU model [5] , so it is not necessarily 
an independent check of the basin structure. Our goal in presenting it was mostly for 
visualization and intuition, since the STU model with its many scalars and charges and 
parameters is rather unweildy. We have not investigated the question of robustness of the 
basin structure in other systems. These include other systems that allow Toda-like solutions 
(cf. [lj's appendices) as well as the large class of Einstein-Maxwell-Dilaton theories that one 
can cook up 12 that are not necessarily integrable in this simple way. 

Recently, there has been some discussion where the attractor mechanism has been con- 
nected to gauged supergravity and also to the possibility of anisotropic geometries [29, 30J. 
It will be interesting to investigate similar questions as we have done here, in these contexts. 

A comment worth making is the connection with supersymmetric attractors. It seems 
plausible that the discussion here can be presented in a form that has analogies with the 



10 See 25J for a discussion of separability of black holes in string theory. 
11 See [22] for an early use of the Harrison transformation. 



12 Even though many of them might not be interesting from the point of a top-down string the- 
ory/supergravity picture. 
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supersymmetric discussion. The solutions here have some moral similarities with the susy 
discussion [311 [32] . 

A question that we have not addressed is that of stability To the extent that non- 
supersymmetric attractors are stable, we believe the discussion in this paper is safe. But it 
will of course be interesting to investigate the stability of these solutions in the full theory. 
It is worthwhile noting here that what we have shown here and in [5] amounts to a certain 
kind of perturbative instability of the (extremal) subtracted geometry. This is indeed what 
it means to say that the subttractor is a boundary between two classes of solutions. 
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A General Subttractors 

Here we present the explicit forms of the metrics and scalars on the various attraction 
boundaries and make some comments about them. To avoid repeating long expressions too 
much, we first define 

H a (r)=Q a (l + d a r). (A.l) 

Note that this is a linear polynomial in r, our solutions will essentially be made out of powers 
of them. In what follows, all non-zero c? a 's are assumed positive. For comparison with other 
papers, we define the warp factor A, as fe 2 (r) = A(r) 1//2 . 

O^T^O case. 

• d\ = d>2 = d 3 = 0. 



Ql Q3 Ad>o Q2 Q3 4flk Ql Q2 2 1 (Q1Q2Q 



3; 



-1/2 r 2 



(A.2) 



In the language of [7] this corresponds to working with an asymptotically linear 
warp factor in the metric (ie., A ~ r) and is the standard (extremal) subtracted 
geometry [5]. 

d 1 = d 2 = 0, d 3 ^ 0. 

4*1 = H s e 40 2 = Q2 H 3 4(fe = Qi Q 2 a 2 = 1 (QiQ 2 )~ 1/2 r 2 
" Q 2 # 4 ' Qi# 4 ' "# 3 #4' "4 ^HJU ' 

The warp factor here goes as ~ r 2 in the asymptotic region. 
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di = 0, d 2 ^ ^ d 3 . 

44* = Ql H 3 40 2 = H ? g 3 403 = 9l H2 n 2 = 1 (Ql)~ 1/2 ^ , 

H 2 H 4 } Q x Hi H 3 H A ' 4 y/H 2 H 3 H A ' 1 

The warp factor here goes as ~ r 3 in the asymptotic region. 



• None of d\, d 2 , d 3 are zero. This is the flat space extremal solution (4.1 4.4), not a 
boundary. For generic positive values of d a these are hairy flat space black holes. 
The warp factor goes as ~ r 4 . 

<) c?4 = case. 

• d\ = d 2 = d 3 = 0. This "vertex of attraction" solution is covered in Appendix [Bj 
The warp factor is a constant, ~ r°. 

• d\ = d 2 = 0, d 3 7^ 0. 

p ^_QiIh r4fe _ Qi H 3 4 , 3 _ Qi Q 2 a _ l (QiQ 2 Q 4 )- 1/2 r 2 

" Q 2 Q 4 ' " Qi Q 4 ' " H 3 Q± 4 ^ • lA ' aj 

This solution also has a warp factor that goes as ~ r like the standard subtracted 
geometry, but the scalars that support it are not of the form adopted in [7, 8] . 

• d 1 = 0, d 2 ^ ^ d 3 . 

40! = Ql H 3 40 2 #2 #3 403 = gl g 2 2 1 (QlQ 4 )~ 1/2 r 2 

"tf 2 Q 4 ' ' QiQS ~H 3 Q A ' a 4 • 1 j 

The warp factor here goes as ~ r 2 . 

• None of di,d 2 , d 3 are zero. 

40x = H i H s 40 2 = H 2 H 3 403 gi ff 2 2 1 (Q 4 )~ 1/2 r 2 

H2 Qa #iQ 4 ' ~H 3 Q 4 ' a A^HJhm: [ ' 

The warp factor of this subttractor geometry goes as ~ r 3 . 

The asymptotic structure of all of these generalized subttractor metrics fall into the 
general form 

„2-ra „ 

» 2 _u - 

o2— n 



ds 2 = dt 2 + -^dr 2 + ff„ r n ^ 2 , (A J 



where cr n is a constant. It is easy to translate this into a more conventional radial coordinate 
R 2 = a n r n , and we find 

n 2 



ds 2 n = —dt 2 + -^dR 2 + R 2 dfl 2 (A.9) 



On 

This is a conical box metric. The cases we care about here are 2n = 1,2,3,4. The case 
2n = 4 corresponds to flat Minkowski space, In = 1 is the conical box of [8]. It is easily 
checked that the geometry has a curvature singularity at R = for all 2n ^ 4. 
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B The AdS 2 x S 2 solution 



When all the d a = in the STU solution we considered, the solution takes the form 

e 40! _ Qi Qz e 40 2 _ Qz Qz e A4> 3 _ Qi Q2 ^ B ^ 

Q2 Q4 Qi Qi Q3 Qa 

a 2 = , = -z. (B.2) 

The explicit form of the metric is then 

ds 2 = ~—,dt 2 + ^dr 2 + R 2 dn 2 (B.3) 

which is precisely of the AdS2 x S 2 form. This is the Freund-Rubin type reduction of our 
original theory where the sphere is supported by the gauge fluxes. 

This solution is therefore qualitatively different from the case where all d a are equal to 
each other but different from zero (and positive). Those solutions are asymptotically flat. It 
is natural to call this solution the vertex of attraction because in the non-negative ^1-^2-^3-^4 
orthant, this solution corresponds to the origin. 



C 4D, Static, Extremal, String Theory Black Holes 

The general (up to hair) four-charge static black hole in M = 4 string theory in 4 
dimensions was constructed in [221 123] . We will follow the conventions of [8] . This solution, 
in a form that is convenient for us is presented in eqns. (A.2-A.7) in [5]. Our charges Q a are 
related to the parametrization used in [81 15] via 

Q. = ^f 2 ^ (C.l) 

which is slightly different from the normalization for the charges used in [3]. The formula 
for the mass is 

M = jJ2 cosh26a - ( c - 2 ) 

a 

We will be interested in the extremal limit, where m — > 0, S a — > 00 while the charges Q a stay 



finite. In terms of the attractor ansatz, this solution is of the form (4.1 4.4) with d a = 2\q~\- 
The gauge field strengths take the form 

F 1 = 2Qi sin 9d6 A d(j>, F 2 = . 2 ®\ — dt A dr, (C.3) 

(r + 2Q 2 ) 2 

F 3 = J= x = 2Q 3 sin 9d9 A #, F 4 ee T 2 = - - dt A dr. (C.4) 

(r + 2Q 4 ) 2 
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D Slicing Through the Attractor Basin 



In this Appendix, we will do a near-horizon perturbation theory of the general solution 



in the single scalar toy model (5.2 5.3). The solution that we perturb around will be de- 
scribed by some (non- negative) parameters d\ = d\ and c?2 = d 2 , and the perturbations that 
are regular at the horizon are described completely by a single 0(r) scalar perturbation 
parameter fa. But since the perturbed solution should also be of the form (5.2 5.3), the di 



and g?2 of this perturbed solution should be expressible in terms of d1,d 2 and fa (See eqns. 



(D.l D.3)). We can determine this relationship by doing a near-horizon expansion of the 
general d\-di solution and matching it with the perturbed (by <f>\) near-horizon equations 
around the d\, d 2 solution. 

In effect, d\, d 2 act as initial conditions for the solution flow in the d\-d2 moduli space, 
and the flow curve is parametrized by fa. Our goal in this section is to determine this 
curve. We will conclude by presenting the full plot of the solution flow along the attraction 
basin corresponding to perturbations around the subttractor solution considered in [5j. This 
latter solution is a specific choice of d°, d 2 , and in figure 2 that we present, the various curves 
are various points in the attraction basin through which the perturbation from the (dj, d°) 
solution moves. 

By demanding that they satisfy the equations of motion, the perturbed fields that are 
regular at the horizon can be expanded in a series as 



a(r) 



a 



d? + 3 d° 2 



+ r 



5 df + 6 d° d° + 21 df 



32 



+ 



y/3 d? fa - y/3 d° 2 



(D.l) 



b{r) 



far) 



d? + 3 d° 

~~ 4 



r 2 / 
32 V 



3 d\ 2 + 6 d{ dl 



3d° 2 



+8^ d? fa - 

-yft d? + y/3 d° 2 



8V3 d° 
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+ 



(D.2) 



+ 



+ 



2 



where 



3 3/8 



a = 



(D.3) 



(D.4) 



2(Q 1 QI)V4- 

If one chooses to work with numerical solutions and start integrating these from the horizon, 
it is useful to have more subleading terms in these expansions for those values of the pertur- 
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I 



-1 



-2 



-3 




Figure 2: The <fi-r curves corresponding to the slice through the attraction basin for our toy 
model when Qi = 7, Q 2 = 4. This is the slice in which the subttractor presented in [5] lies. 
As we have emphasized, the qualitative features of the plot are not different from our rather 
ad-hoc circular slicing of the basin presented in the main body of the paper. 



bation <pi for which the 0(r) term in the scalar expansion vanishes. It is straightforward to 
determine them by demanding that the equations of motion are satisfied order by order in 
r and we have computed them, but they are too cumbersome and we will not present them. 
In any event, since we have the full analytic solutions, we can always translate any question 
we want in this context into an analytically tractable question. 

The near-horizon expansion of a general solution with parameters d\, d 2 is 



a(r) 



b(r) 



di + 3 d 2 
4 



+ r" 



5 c^ 2 + 6 d 1 d 2 + 21 d 2 z 



32 



1 

a 



1 + r 



d x + 3 d 2 



2 

+ h ( _ 3 di * + 6 di d2 



3d 2 2 



+ 



+ 



(j>{r) 



VS. / Q 2 



+ r 



-\/3 d\ + %/3 d 2 
4 



+ L ( 3 ^3 d x 2 - 3^ d 2 2 ) + 
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(D.5) 
(D.6) 
, (D.7) 



and we want to match this with the previous perturbation expansion to determine d± and 
d 2 in terms of d®, d 2 and <f>\. It is easy to see that the entire series expansion matches term 
by term when we set 



d\ = d\ — VS(f>x, d 2 = d 2 + —j= 



(D.8) 
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Eliminating 0i (which is essentially the parametrization of the curve in the attractor basin), 
we get the slicing curve (in fact a line) through the d\-di attraction basin to which a black 
hole with parameters d\ and belongs to: 

d 2 = \(d° 1 -d 1 ) + d° 2 (D.9) 

As final illustration, we plot the flow in the attraction basin (in 0-r space) starting with 
the subttractor solution we considered in [5], which corresponds to 

We can do this by numerically integrating the perturbations carefully as we mentioned before, 



or by plotting the various solutions along the curve (D.9) with d\ and d® chosen as above. 
The result is plotted below and as we emphasized, the qualitative features are identical to 
what we saw in the main body of the paper. 

Similar constructions can be done for the 3-scalar solutions in the STU model as well. We 
won't present the details of the perturbation theory because it is merely more complicated 
and not conceptually different. But it turns out that the hyperplane that cuts through the 
attraction basin in 4-4-4-4 space under perturbations of a black hole specified by d®, d®, 
g?3, g?4 is given by the simple relation 

d 1 + d 2 + d 3 + d 4 = d^ + d^ + dl + dl (D.ll) 

Note that in the STU model there are three independent scalar perturbation 1; 2 , 03 now, 
which correspond to different possible directions on the hyperplane for the solutions to move. 
To obtain the equation above, we have eliminated these perturbations from 

C?l = d? + 01-02 + 03, 4 = ^- 01 +02 + 03, (D.12) 
4 = ^3 + 01 + 02 - 03, 4 = ^4 - 01 - 02 -03- (D.13) 



These are the analogues (D.8) in the STU system. 
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